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ABSTRACT 



The state variable formulation of a Foster Wheeler 
ESD-fflboiler is developed from fundamental principles. 
The' response of the model for various input signals is 
determined using CSMP -I, the IBM simulation language. 
The sensitivity of the model to various coefficient 
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NOMENCLATURE 



As , Ay, A B = superheater, riser and downcomer cross-sectional 



z 

areas respectively (ft ) 

C = evaporation level proportionality constant 
(ft-sec/lb) 

C S >C,,C„ = heat capacitance for superheater, riser and downcomer 
tubes respectively (BTU/lb-R) 

C a = average heat capacitance of combustion gases times 
air- fuel ratio (BTU/lb-R) 

C; = heat capacitance for feedwater (BTU/lb-R) 

C^ s = average heat capacitance of superheated steam 
(BTU/lb-R) 






o 



3 

h 



6 “ 



friction coefficients for superheater, riser and 

2 5 

downcomer tubes, respectively (sec /ft ) 

2 

acceleration due to gravity (ft/sec ) 
enthalpy of saturated vapour corresponding to P 

B 



(BlU/lb) 

W8= enthalpy of saturated liquid corresponding to P 

B 



( BTU/lb) 

H wt) = enthalpy of liquid in mud-drum (BTU/lb) 

Kn = enthalpy of drum and downcomer liquid (BTU/lb) 

enthalpy of evaporation corresponding to P (BTU/lb) 



K 



95 ' 







air cooling coefficient at superheater bank (BTU/sec) 
air cooling coefficient at riser bank (BTU/sec) 
heat-transfer coefficients from combustion gas to 
superheater tubes, and from superheater tubes 
to steam, respectively (BTU/lb-R) 

heat-transfer coefficients from combustion gas to 
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K e = 



i K & = 



**S) L V>**D — 



M = 
Ms = 
Mb = 

P* = 
Ps = 
P w = 
a s = 

Ss = 

= 

^ 90 = 

T( = 
Tb = 

T S = 
T W = 
T$w = 
T aw = 

^ 3 a- 

T 9 b = 
T c = 



V 



riser tubes and from riser tubes to boiling 

li'guid respectively (BTU/lb-R) , (BTU/lb-F ) , (BTU/lb-R) 
evaporation rate constant of drum liquid (lb/sec-R) 
constants for state equations of saturated steam 

o 2 -» 

( R-ft /lb) , (ft ) 

superheater, riser and downcomer tube lengths 

respectively (ft) 

mass of drum liquid (lb) 

mass of superheater tubes (lb) 

mass of riser tubes (lb) 

2 

drum pressure (lb/ft ) 

2 

superheater outlet pressure (lb/ft ) 

2 

mud-drum pressure (lb/ft ) 

heat-input rate from tube walls into the superheated 
steam (BTU/sec) 

heat-input rate from hot gasses into superheater tube 
walls (BTU/sec) 

heat-input rate from riser tube walls into boiling 
liguid (BTU/sec) 

heat-input rate from hot gasses into riser tube walls 
(BTU/sec) 

o 

feedwater temperature (R) 

9 

saturation temperature corresponding to P (R) 

B 

o 

superheater outlet temperature (R) 
drum and downcomer liquid temperature (°R) 
superheater tube-wall temperature (°R) 
riser tube-wall temperature (°R) 

average gas temperature at superheater banks (R) 
average gas temperature at riser banks (R) 
comnustion gas temperature entering superheater banks 

o 

<R) 

3 

volume of vapor phase in drum (ft ) 
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V y = velocity of riser mixture (ft/sec) 

V w = velocity of downcomer water (ft/sec) 

3 

V = total drum volume (ft ) 

W s = steam mass-flow rate at the superheater outlet 
(lb/sec) 

Wf = fuel mass-flow rate (lb/sec) 

Wi. = feedwater mass- flow rate (lb/sec) 

W w = downcomer mass-flow rate (ln/sec) 

W = riser mass- flow rate (lb/sec) 

W a = air mass-flow rate from blower (lb/sec) 

W ac = chemically correct + 50 % excess air rate (lb/sec) 

W e = mass-evaporation rate from drum liquid surface 
(lb/sec) 

W B = steam mass-flow rate from drum into superheater 
(lb/sec) 

X = quality of mixture leaving riser 
Xy = throttle opening (%) 
y = drum liquid level (ft) 

3 

/> B = saturated vapor density corresponding to P (lb/ft ) 

B 

3 

y°s = superheater outlet density (lb/ft ) 

3 

>°w = saturated liquid density corresponding to P (lb/ft ) 

B 

/> = density of liquid vapor mixture leaving riser 
(lb/ft 3 ) 

)) = evaporation level (ft) 
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I 



INTRODUCTION 



Boilers as understood by marine engineers are closed 
vessels containing water which by the application of heat is 
converted into steam at any designed pressure. This steam 
is than used for the production, through machinery, of 
useful work. A dynamic model of steam turbine machinery 
consists of a boiler model and a turbine model. Tne 
difficult part is the boiler, as a load change causes 
variations in some important properties such as toiler 
pressure, temperature and drum water level. The analysis in 
this paper will use Chiens dynamic analysis (1) of a boiler 
as a ref erence .Chi en considered a naval boiler which for 
purposes of analysis was divided into four sections namely a 
superheater, a downcomer- riser loop, a drum and a gas path. 

The principles of thermodyna raics , heat-transf er and fluid 
mechanics were used to describe the dynamic behaviour 
corresponding to each section of the boiler and these were 
derived from equations of continuity, energy, 
heat-transfer, and momentum. The equations involve partial 
differentials as well as nonlinear terms. These equations 
were reduced to the ordinary linear equation form by 
applying small perturbation and difference equation 
techniques. Linear equations thus obtained were reduced to 
ten state variable equations and solved by digital computer 
techniques. 

Since there is an increasing interest in boiler 
modelling, the objective of this thesis was to develop a 
comprehensive boiler simulation model in a form useful for 
modern control (i.e., multivariable control) analysis. 
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II. 



BOILER MODEL 



A. BOILER CONSIDERATIONS 

The control problems of high pressure boilers have 
become more and more critical, both from the operational and 
the economical points of view. A dynamic analysis is the 
method to be used for a control- system analysis. It 
consists of a complete understanding of the process to be 
controlled and the effects of physical and chemical changes. 
The analysis is not exact by any means but the results 
obtained should be in good qualitative agreement with actual 
tests. The major difficulty in boiler analysis is the fact 
that the whole system is very complex and contains numerous 
variables. 

Chien(l) considered a naval boiler which for purposes of 
analysis was divided into four distinct sections namely a 
superheater, a downcomer-riser loop, a drum and a gas path. 
The detail of analysis is described in Part B of this 
section. The basic equations used in the analysis are those 
of continuity, ener gy (heat-transfer) , momentum and the state 
equations. These equations involve partial differentiation 
as well as nonlinearities. Generally , the equations have the 
form 

j (X , y , 2 , ) = 0 

\ 

To eliminate the nonlinearities one uses perturbation 
theory, which effectively approximates the response of the 
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system to small signal changes about a chosen operating 
condition. Thus, the equation is perturbed about its steady 
state operating condition to give the linearized form. 

Hence it can be written as 

^ AX +• Ay + -&L A z + = 0 

ox oy 62 

The perturbed variables are AX, AY, AZ, etc. and the 

partial differentials ( d£, 2>f, df, etc.) that form the 

ax liy Th 

coefficients of the perturbed variables are evaluated at 
steady state operating conditions. This technique was also 
followed by Whalley (2) in modeling the same boiler plant. 



1 
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Figure 1 - A CROSS SECTIONAL ARRANGEMENT OF THE FOSTER 
WHEELER D-TYPE MARINE BOILER 
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B 



BOILER ANALYSIS 



The boiler studied in this analysis is a Foster Wheeler 

D-type- marine boiler; which is an oil fired, two-drum, natural 

circulation unit having a rated output of 28,800 lbs/hour at 
2 

350 lbs/in gauge, 1200 F. A cross sectional arrangement of 
the unit is shown in Figure I.The following assumptions 
apply to the physical simplifications in each of four 
sections. 

1 • Superhea ter 

(a) The inertial effects of superheated steam are 
neglected . 

(b) The superheater tubes are assumed to be a single 
capacitance w-ith restriction on the drum side and another 
restriction on the load side. 

(cj Desuperheaters are not considered. 



2. Downc ome r rise r loop 

/ 

(a) Only natural circulation exists. 

(b) No boiling takes place in the downcomers. 

(c) Vapor and liquid velocities in the riser are 
identical. 

(d) Heat-transfer rates to the boiling liquid from the 
tube wails are proportional to the cuoe of temperature 
difference between the wall and the liquid. 

(e) Steam quality is uniform in the riser. 

(f) Liquid temperature is always the same as the 
saturation temperature corresponding to drum pressure. 

(g) Downcomer liquid temperature is the same as the 
drum liquid temperature. 
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3 . Dru m 



(a) There is no temperature gradient across the drum 
vapor phase^and the temperature is always the saturation 
temperature corresponding to drum pressure. 

(b) The liguid phase has no temperature gradient other 
than across a very thin boundary layer at the drum surface. 

(c) Evaporation or condensation rate in the drum is 
proportional to the difference of liguid and saturation 
temperatures. 

(d) Feedwater temperature is assumed to be constant. 

(e) Liquid-level changes due to bubble formation in the 
drum are neglected. 

4 . Gas Path 

(a) The air -fuel ratio is assumed to be constant. 

(b) Temperature of combust ion gas entering superheater 
is proportional to the firing rate. 

(c) Waterwalls are lumped with the riser-banks. 

(d) The heat-tr ansf er rate at each tube bank is 
determined by the tube wall temperature and the average gas 
temperature. 

(e) Inertia of the hot gases is neglected, that is, 
velocity changes take place instantaneously. 

(f) Delays due to the heat capacitance of the hot gases 
are neglected, that is, temperature changes take place 
instantaneously in combustion gases. 

(g) All heat transfer is due to turbulent convection 
and radiation. 

The following steps are taken in developing the 
equations for each of the four sections. A simple schematic 
diagram of each sub-section is included at the beginning of 
each part of the analysis. A brief statement of the physical 
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situation is included under the headings. The dummy 

coefficients such as a, b f C , D , were used for 

1 1 11 n 

convenience^and their values are declared in the Appendices. 

the quantities such as W , h , "h , etc. are the steady state 

B S B 

values of W , h , h , respectively. These values can be 
B S B 

evaluated from the original unperturbed nonlinear equations 

by setting all derivative terras to zero and solving for the 
unknown values in term of known values. 
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Figure 2 - SCHEMATIC DIAGRAM OF SUPERHEATER 
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5 



Superhea ter conti nuity 



Mass flow in - Mass flow out = Hate of change of enclosed 
mass 

w„-w 3 = ic-sv.) 

at 

The perturbed eguation is : 

A w 6 - A W $ = L s A s SA/> s 

or, 

AW B + a aw 6 > O a SA/= s (i) 

where S is the [ Laplace variable denoting differentiation 
with respect to time 

6 . S uperhea ter energy b ala n ee 

Heat input from tube wall + Heat in incoming steam - Heat 
with outgoing steam = Rate of change of internal energy 

Q+wk-wh = i_(L A bk ) 

8 Be > s s ^ ' b °J s s; 

The perturbed eguation is ; 

- L s\(/ 6 SAk s +R s SAy> 5 ) 

Since,Ah &3 0, ^k^C^ATg and L s A,.SA/ s = AW B ~ A W § 

After substitution and rearrangement: 

( L s A sA c fv*) SAT s+ < fi s -k B )AW B - Aa s --w s c„,ATs (2) 

or. 



a 4 SAT s + a s AW a + a 6 AQ s - a 7 AT S 
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7 . 



H eat conduct ion across s uperheater walls 



Heat input to wall - Heat output from wall = Rate of change of 
energy stored in superheater wall 



Q 



as 




The perturbed equation is 




A&gs ~ Aft s = M $ C s SAT sw 



or. 



A ®-<js 



a Aa 

9 S 



= a s 

io 



At, 



Sw 



( 3 ) 



8. Heat input to supe rh eater wa ll 

The following empirical formula includes the cooling effect 
if air supply W^differs from ideal air supply (whalley (2)): 

Q aS = K 95 W r(V- T S »)-Kas(l-V ) 2 

Wac 

The perturbed equation is 



A a = 
9s 



( Ts5 - T s«) 4 ^ * + >4W 

Wac 



or. 



a 40 

n 



9s 



a 12 Avjj + a 13 (A T gs - AT SW ) + a 4 Aw a 



( 4 ) 
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9. Hea t tran sfer acr oss s uper heater wall 
For turbulent gas flows: 

a s ■ W Vw-Tj) 

The perturbed equation is : 

Aa s ’ ^T. 2 K S^SW-T 5 )AW b + K 3 W b 0 ’ S ( At&w _ At j 

or. 



a i5 Aa s - a i6 AW * + a l7 (AT sw - AT S ) 



C 5 ) 



10. Heat tr ansf er from com b ust i on gas 

Average gas temperature at superheater wall is T where 

3s 



9S 



T c - 



0-5 

s w f 



and is related to the heat capacitance of the combustion 

gas by ; 

C 3 

where C is constant 
c 

Thus ; 



- 2C 0 (w f+ Vl a )( Tc -T ) 



The perturbed equation is ; 



Rearranging ; 

^ i& , 5 - ( 5 'f^\)iT c + CT e .T 3s )M f + (f £ -f 9s )AV(Wf^a)4T 9s 

or. 



“i* L S* “ a ia iT = + a » iw t + a 21 AW « + W 



(«) 



1 1 . Superheate r momentum 
Neglecting the inertia of the steam ; 



( p *- p s K-(^)v B -« 



Since : W 



& 



*A A sV e 



<v 



= 0 



Substitute : 1 f L s 

5 * 9 \kl 

The perturbed equation is ; 



AP ~ A P S = 2 ls — Aw a - i- — Ap 
& S S B 7 * A? /* 



A 



Rearranging ; 



or. 



2i s |- a Av( B +AP s -f s ^ 2 Ay> & = a P, 

V B A 



& 



Q 23 AW & + a 24 AP 3 +a 2 5 A A = a 26 AP & 



(7) 
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w 



w ) 



A 



Figure 3 - SCHEMATIC DIAGRAM OF RISER-DOWNCOMER LOOP 
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12. Riser con t inuity eq uati on 

Mass flow into riser - Mass flow out = Hate of change of 
enclosed mass 



A ( 

W - W . *- 
w dt 

The perturbed equation is ; 






Aw u - AW = 1 



A^Ly.SA^P- +A y ^sA/>] 



since 



then. 



± = _X + Cl-X) 

7° A A 






If the starting point of evaporation in the riser tube is \j 
above the water drum exit where ; 



Y) * 1 - CW w 



then. 



A>) - - C AW w 

A v ^ S Ay° can be neglected since ^ is small 

Substituting ^ and A ^ and rearranging ; 



-2 



0 • 5 - -4- |fiA»4.0,S A..I..<K I 



SAX + 0.5 AyLyX K B [ ]sAP ft 



+ 0-5 (y5 + f>jCS A W w 



= Avj w - Aw 



or. 



a 7 s AX + a 2QS AP & t a 29 SAW w 



— ^30 * ^3i AW 



( 8 ) 
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13. Riser momentum equ ation 



Pressure drop - frictional losses - gravity head - 
exit loss - monentum flux change = Inertia force 



YJ 



-•P A v - 



2 

Ww 



1CAAv) ; 



= JLy dw 

1 dt 



Therefore : 



_A^q( Dv> + 2 



1 



L C j 

W Ww _ l /} L v* d w 

/ 1/wA* v7 ’ 



After perturbation and substitution for AyO as done before 
in the Riser Continuity Equation, the perturbed equation is 



(ap„-ApJ-2 F [(4- £)5*+£]aw-f»*(1 



iW ' FW '\'A) AX 



_ _2 

+ T™. - L 

Ay. 

where 



JL£_ + o 

V 



F = 



AVb 

( 



(4-4-i 

V/>w y°a J 



AX = 



fv W + _3 



Bearranging ; 






waere 



ap w - ap b - (\ + a v s) aw+ b x Ax- 

b * * 25 (^ + f)/f A 7 3 



L v 



f 



SAW 



2W„ 



f a;a w 



AW + c< 



% 



a v = - 2 L 
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and is small enough 
to be neglected 




oC = 






Therefore ; 



or. 



4P,-AP t = k AW + a SAW + b,AX- 1*! Aw, 

■J 



a 3 ,AR» + a „AP 6 = o 34 Avi+ a J5 SAW + a 36 4x + a 37 AW (9) 

14. Hear balanc e acro ss riser tube wall 

Heat from hot gas to riser tunes-Keat from tubes to mixture 
=Rate of change of internal energy stored in tubes 

o 3 ,-“* -df (H * c ‘ TM 

The perturbed equation is *, 

' Aa e ‘ M 6 C 6 SAT 6( , 

or, 

a A a QIS + o A a -a. sat 

38 39 B 40 U feW 



( 10 ) 
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1 5 . Riser heat balance equation 

Heat input to risers- Heat output from risers 
=Rate of change of internal energy 

V ,I A»- Wh • | t [o-«A f (i-*-i))(/>y>,)k] 

From the continuity aquation ; 

j. [°- 5 = w w -w 

After substitution and perturbation ; 

A \ + \ All w I , + ^v J t» _ ^)AW w -W Ab = 0.5 (/+/>,) S Ah 

since : 

k = k Wft +xk n 

If changes in latent heat Ab^are neglected , then 

Ak=Ak WB + b^AX 

From steam tables ; 

A 11 -! > 

A b = At w 

Therefore ; 

0.5 A,L„(^+^ w )k c SAP b + 0-5 SAX - Aa B 

= (b WI) ' K) A\ - wk c ap b -wTi^ ax + w w at w 



or. 



a 4l SAP ^ 



+ a SAX+a, Aa * 

42 43 B 



a 44 AV 'w + a 43 AP B + a 46 A X +Cl + 7 AT W 



(il) 
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16. Riser heat trans fer equatio n 

An empirical formula that accounts for turbulent gas flow 
is adopted and gives J 

=V W f (V~ 'bw)“ + k v( t ^b- t m ) 

The perturoed equation is ; 

Aa gB = Aw j + Kj.Sf (at,.-at.J 

w f 

+ ^ AVJ * f AT,.,) 



_ 2 

Vi 

ac 



BW 



Rearranging ; 



i4 3»' K 9» W * AT 9» - ' K 3 si f at m + — AWj 



Aw a ♦ 4 K rT ,*. AT, 

Vi a.c 






- 4K r T 0 „ AT bw 



hence ; 

‘■33 V ’ gE, 



Aa <i»'K »^* +4 K -V) AT r 



0-6 



--(K,»W( +4K,f s „) AT„ 

(.TgB ~Trw) 



— 0 . 4 . 

Wf 



■A^j. 



4. ^ Ka.y* ( W ac - Vi o. ) ^ 

A ^ 



— £ 

w 

" ac 



Of/ 



a 48 + a 49 AT ^ = a Sl3 AT B w + a 51 ^ w j: + a K , A W a 



(12) 
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17. Rise r heat trans fer to bolii nq fl uid 

The heat transfer rate is assumed to be proportional to the 
cube of temperature difference between wall and mixture. 
Continuing to follow Whalley(2): 

Q & = ^ T bw ~ T r ) 

The perturbed equation is ; 

h \ = 3yT M -f,HAT B „-A Te> ) 

cr, 



a S3 ^ a e. “ ^S4 at bw 4 a ss 



(13) 



18. D owncomer m omen tum equat ion 

22 Forces due to pressure difference-frictional loss 
+gr avitat ional head- entrance loss = Inertia force 
Hence ; 

(P B -r,K_ ( Vw)yj, At, A p (Lp+2) — 

u h ^ u 

d (V»AM 
dt <j 



since * 



-A a j> v w 



<1 o 

w w , ^ , Ww 



(p B -p„)- ~ •" i — 

V» 6 .A» 0 



(p b -p w )- fitb> + o.5) _L_ w»f/>L _ h i.(v»w) 
^ R ; /gAAfJ * gA 0 at 



= d _U Wv 

2 $Vw dt 

2 • cl 
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The perturbed equation is ; 



A A P w * ( b w + Q w S ) A W w — L d Ay° N 



where : 




Since flow is incompressible ( y0, v is constant both in time 
and space), then - o 

Therefore ; 



19. Drum steam mass balance eq uation 

Y Vapor rates to/from drum = Rate of change of vapor mass 
in drum 
Thus } 

Evaporation rate + vapor rate from risers-steam 

rate from drum = — ( V a P a ) 

dt 7 

Therefore : 




(14-) 




The perturbed equation is ; 



A^e + X Aw + w AX - Aw fl = V 6 S + /\S A V & 



Substituting for ; 



A V a = - A Ay 



A A a K ftAP 6 
Aw e - k e (AT w -rr B ) 
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gives ; 



V B K B SAP e A S A y « K e ( AT w ~ AT 6 ) + X AW + W AX - A W fi 



oc. 



a 6c^AP B + a gi S Ay = a 62 (at w -at 8 ) f a 63 Aw + a 64 AX + a aw s (is) 

20. Dru m l iq u id, heat balance equation 

£ Heat to and from drum with liquid = Rate of change of the 
internal energy of the liquid 
Hence ; 

Internal energy of riser liquid + heat transfered in 
feedwater - heat taken in by downcomer liquid - latent heat 
released with evaporating liquid =_cl_(drum liquid internal 
energy) ut 

Therefore ; 

C v u- x) WT S + C; W;Tj - - W, h BW = i(C D MT„) 

at 

But C y ) C; ^ Cj 3 i ■ 

The perturbed equation is ; 

( i~ X )T fe A W + ( X ) W AT B - WT B AX + Wt AT ( ^ A W £ - W w A T W 

-Ty, AVJ w -Vl e AK 0w -.h BW AW e = MSAT w +T w SAM 

Substituting for ; 

AW e = K c ( AT W - At 6 ) 

^ =• K c A P & 

AM = Ay 

and T. f* constant, so AT. ss 0 
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Thus: 



(l-X)T b AW+ (l-X)WAT 8 - WT a AX +Tj AW— W w AT W 



”* AWfl - W C K C AP a 



-^ B „K e (AT w - AT 6 ) 



M S A l w + A^ W T W SAy 



or. 



a 66 aw+ a 67 AT B + a &s AX + a 6g AVJ t f a 70 AT w +a 7i AW w + a 72 Ap B 
^ ^73 ( ^ 'w “ ^ t b) — *^74 S AT W +- a. 7 _S Ay ( 16 ) 



21. Dr um l iq uid mass balan ce eg ua tion 



Liquid input/output rates 
of liquid 
Hence : 

feedwater input rate + riser 
liquid rate-evaporation rate 

Therefore • 

# 



= Rate of change enclosed mass 



liquid input rate-downcomer 

= ( M ) 

dt 



w j + ( 1- x ) w - W w - w c 

The perturbed equation is ; 




AW; + (1- X) AW - W AX - AW w - AW e s S AM 



Substituting for : 

AW e = ( AT W - AT B ) 



AM a Ay^yy A y 

gives : 

-K e AT R + A^P w SAy = (i-X)AW-wAX-AW w -l<eAT w + AW t 



a 66 ^B + a 7, s * y 
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22. Equation from throt tle valve 

'"s “ h (*».> p » ) T J 

The perturbed equation is : 

AW = A X v 4- A P s + AT 
s Ax Ap s s At s s 

or, 

°9.1 A W S - fl 94 AX v f a 9s 4 P s + a 9(j i T t 

23. State e quati ons 

. At b . |l? AP a - K t ap s 

or, 

a t 6 - a 97 ap b 

B = ^ P 6 = K B ^ P & 

or, 

A/a = Ss A P & 



Cl 8 ) 



C i 3 ) 



C 20) 



A _ A T ga 

Akw=-^“AP a _K c AP 3 
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( 22 ) 



or, 



A h 



BN « CX g g A P B 



p* 

A P S 



‘ f*(Ws) 

.^ iTs + |iSA A 



or, 



A Pi = a io^ T s + a i°i A/s 



(23 ) 



or, 



A T c - Kf AWf 



AT r - a L0 2 



( 24 ) 
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C. STATE VARIABLE EQUATIONS 



The following procedure is an algebraic method to reduce 

the twenty-four governing equations to ten state variable 

0 

equations, that is, in standard matrix form of X(t) = AX(t) + 
BU(t) 

Here X (t) is the matrix of system state variables and 
U (t) is the matrix of system inputs. 

From (1) 



a A 



s /« 



4 A 



= a, AW -i- a AVL 

15 2 => 



- ib AW Q 4- — Aw 



a ; 



6 a 3 s 



or. 



A A = b i AW b + b 2 A 



W 



( i.l) 



From (2) 



From (5) 



a + AT s 



a 3 “ a e Aa s + a 7 At s 



a j s .isAw At. 

s 0 . 4 . s a 4 s a 4 ^ 



(i-2) 



Aft . Aw a + 2H 4T a _ Si iT s 

s aio 0-13 



(1.3) 



Substituting Q from (1.3) in (1.2) and rearranging ; 

s 



At. « 



| a l6- a t7 ! 


0^ 


1 AT c- 


( 

— + 


Alt- Ql6 


\ a 4-*i8 


^4 J 


I $ 


l°* 


a i ' a i5 j 



AW - — ■ At 



sw 
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(1.4) 



or. 



AT S - V T s - b 4 ^ 8 -b.AT 



* $w 



From (3) 



o»s4t„ . a a j4 + a 3 a a s 



From (4) 



At . 3 ii Aa + £i AQ 
s " a i0 <js ai0 s 




^/ w f + a7:V£l iT ~ +a “ iw * 



From 




^19 AT + AWj + ^i- 1 AW a + 

ciia c d-i* fl id 



°-2Z 



Setting (1.6) = (1.7) and rearranging'. 




(1*5) 



C i • 6 ) 



(1-7) 



(i. 8) 



Substituting AT = a AWf 

i°z r 

rearranging : 

la 



AT _ 

3 s " 



13 • ^102 _ a i Z ^ Oio 




from 



O-t* tt u a i« J 



Aw 



i+ 



(24) 



into 



( 1 . 8 ) 



/« 13_«£S 
\^u j 



and 
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or. 



A1 V = b 6 iW f + ^ AT s»- b . AW a 



Cl. 9) 



Substituting ATcj 3 from (1.9) into (1.6) and rearranging : 



Aa = 
9 s 



K + + (k®ii-?>a) 4 T +fa. ^*12 W 

a u a a j ' \ a Li a u / ^ 14 a tl ) * 



or. 



AQ 3S = b 9 Aw ) + b io AT s « + b i, Aw c 



( 1 - 10 ) 



Substituting A Gl , from (1.10) and AQ.. from (1.3) into 
(1.5) and rearranging ; 



At.. -(^2-1 A«j + 



Svg 



a io 



/ q 8- b io t a 9 • Q i?\ AT t / V M Aw 






io ^*10 ‘ ^ 



15 



SVf 



4o 



a 



or * 



+( Vdii \ Aw - ( a > a » ^ at. 



At s , » b x. Av, f *• b i3 lT « * b «‘V K, AW„ - b M At s 

From (17) 



(1- n) 



Ay , 



_ At r ^ Aw -h AX 

0-77 



a 



77 



a 



77 



+. ^ AW, AT, + AW; 

a 77 <*77 <* 77 



( 1 . 12 ) 
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Substituting AT R from (19) into (1.12) 



s 6 . 


*9 7 


ap r + 


a 38 


Aw 


a-77 


o 


fl-77 




Ag c 


AW, 




AT U 


+ 


a 77 




a 77 





a 



a 77 



77 



or. 



Ay = -b 7 AP 6 + b i8 AW + b ig Ax + b 20 Aw w +. b 21 AT w + b 22 Aw- 

From (16) 

AT W = 



a?3 A y + Aw + AT ft + AX 

a 7+ a 7 + a-7+ a7 + 



+ 



AW; + (Aslfzi) AT„ + £» AW w + ^ AP, 
a 74 a 74 a 74- a 7+ 

* 



Substituting Ay from (1.13) and AT r from (19) into ( 
and rearranging • 



A T w _ (* 75 ' bl7 + Q 9 7 • a 67 _ a 97- q 7.9 ] A P R + 



^ _ Q 75- bia | ^ 



a 



74 



a 



74 



a 



74 



a 74 



a 74 



l 6S 



a 75 • bjs | A X + ( — — - 0,73 ' ^ 2C 



a 



74 



K„+ a : 



a 7 4 

^73 • b-4t 



>,-^73 itii 



a 74 . a 74- 

® 75 ■ b,7 ^ ^ 



a 74 






a 



74 



a 



74 



Aw„ 



Aw- 



or. 



A ^ w ' ^3 Ap 6 + b 24 Aw + b 2s Ax + b £ 6 Aw w 4 b 7 AT w + b^AV^ t 



(1.13) 



(1.14) 

.14) 



Ci.15) 
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From (15) 



A p Ay+^ AT W _ iii£ A T +-^AW+iiii Ax+^iiAW (1.16) 

ft a 60 ' cuo a 60 a eo ci&o a 6o B 

0 

Substituting Ay from (1.13) into (1.16) and rearranging *. 



AP 0 




^*U»+(a»- a -^) ix 

a 60 J [ 0-bo cuo ) 



At w 



or, 



AP 6 * b 29 AP B + b 3c AW + b 31 Ax - b , 2 Aw^ 



+ b 3.-> AT w " ^AWi t- b 35 A W £ 



(i-i 7 ) 



From (11) 



» a Al 6 
AX - - A P R 



a 



42 



- Aa + Ali Aw w t ^ - A p £ 






a 4ii 



CU2 



4-^A. 6 Ax + A Tv, 

0 . 4.2 ^42 



(i-is) 



From (13) 



a a = _ AT 



^55 



BW 



AT, 



1 S3 



a» 



E> 



5-3 



40 



Substituting AT B from (19) 



A Q. s At + Qa5 ’ a 

W B rt 6W + - 



a S3 



•9, A p 

u 53 6 



( 1 • 1 9 ) 



Substituting A Pp, from (1.17) and AQ.„ from (1.19) into 

r> 5 

(1.18) and rearranging ; 



AX ^ [ Q l 5 _ q »3- a 53-^97 ^ A p _ 



, a 42 



a 42 



a 



42' tt 5 3 



B 



a 42 



AW 



“«_ j^A») 4x + /£m + iitbjj. 

a 4i ^42 / l a 42 a 42 



A w w + f— -iliii]AT 

\ a 42. <*42 ) W 



^41' 



a 



A W r - 



/* a +l * ^35 \ 



42 



a 



42 



AW 6 _ 



AT, 

q 42 ■ 0-53 



BW 



or, 

AX = ^36 A P B - k 37 AW + b 38 AX + b 33 AW W 

+ ^ AT W +■ b 4i AW( - A W R - At Bvv (t.2o) 

From (8) 

AW , - °2L AX . ii« AP. +122 AW„ + tL‘ AW (121) 

^29 &29 A £g ^29 

e * 

Substituting AX from (1.20) and A Pfi, from (1.17) into 
(1.21) and rearranging ; 



41 



AW„ 



I , \ 

^ Z7 ‘ ^ 






a 29 



Ap a + 



^27' ^37 
a 29 



a 27- ^38 + 029^31, + 

<3-23 & 23 

' / 

A \ + 



^2rt * b 



V 



a* u 32 



29 



Biliks. + 2i‘ | AW 

Q. £9 Cljjg 

^30 ^27- b 



l 27 - 033 



&29 



fl 27 • b 4<3 ®28 • b 33 



a 






29 



< 3.29 



^ 29 • b .* 4 fl 27 • k 






•29 



27 • °41 
a 29 



a 29 '^ Ww 



A W( 



^27-^42 ®2S ,b 33 



a 



2.3 



<323 



A V 



(, -i> 43 \ 



27 ’ 

a 



29 



At, 



BW 



or. 



AW w = -b 44 AP 5 + b 4 - AW - b 4& AX + b 47 AW w - b 43 AT W 



b 49 AW { + b so AW b + b sl AT bw 



(l-22) 



Substituting from (1.22) into (14) and rearranging ; 



ap w Jlmh 3 ‘ ^.) AW / Q ^ a ^/. b JdU P + 



d 59 , b 43 



k 57 






O. 5 7 



a 



AW 



57 



°59 . b 4 6 j ^ x — 



f a 59 ■ 



■57 



0-67 



AT W + 



a 59. *>49 



a 



AW: 



57 



a 53* bso 



a 



57 



aw b + 



^59 , bs 



sa ■ 051 ) AT 

0. S7 



B'rt 



42 



( 1 - 23 ) 



or. 



AP„ = b 32 AW w - b S3 L P R + b^AW-^.AX-bsefiTw 



* b s ,AW. + k S9 AW 6+ b 39 AT 



BW 



From (9) 

AW s !i 2 A P„ + A p* - h± AW _ AX _ £22 AW 

w ^ B rt a as 



a. 



a. 



'33 -33 a a« “3s 

Substituting A P* from (1.23) into (1.24) and rearranging ; 



AW s 



®3; • k s2 U J? 



a 3S Q 3.>? 



- — Aw n + 



°33 a 3’< b 



S3 



a 



35 



■33 y 



AP S + 



B— j & ~ * \ 

— — - r ^- 4 AW 

fl 3o ^35 / 



^3^ ' ^S3 Q 36 

a 



a„ 

3a 



3.5 



\ 






f \ 


AX ._ 

/ 


fl 32 • b 5s 


‘V 


fl 3 2. bs 7 


u a ~ 

3o ) 


a 3S / 



AW; 



«3;.b 33 



AW. 



a 
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^ q 3r • b 33 

a 33 



AT r , 



B'W 



or, 

AW * b 6o A'W w + b 6i AP ft +■ b - 2 Aw - b 6 ^ Ax - b^ at w 

+ b 65 AW. + b 66 AW &+ b &7 AT a 
From (10) 



e>w 



( i- 2 s) 



0 

AT 



&w 



= -il Aa +, ^21 AQ 

CI 40 ( 1 ® a 4o 1 



(1.26) 



43 



From (12) 



a 4* ia 9 6 * a S » iT S »* %1 A V Q S= 4 V- a 43 4T 90 ( L27 1 



Since AT *AT fl from (21) , then substituting AT 

C)5 ^ gs 

From (1.9) into (1.27) and rearranging ; 



Aa 



3 ft 0 . 4 $ 



S0 AT bw + 



^51 ~ & 4-3 • b fo ^ AWr . ( ^sz + *-* 4,9 ‘ ^ 



a 4-s 



A Wi 



a 48 



AW„ 






a 4« 



At 



sw 



( 1-28) 



Substituting A from (1.28) and AQ. R from (1.19) into 
(1.26) and rearranging ; 



AT.. - 4 At _ 4 ( 12 . \ fa sl -a,b) iw 






+ (^/( a “ t ‘ wb v i "« 



. Cl 49 - b 7 

^ 4*0 * 3 - 4 $ 



AT 



sw 



or # 



+ f a 3S- a BS ‘ Q 97 j A P p 

I a 4 o • ^-55 ) 



^ T bw = ^68 AT aw * ^69 A Wj. + b 70 Aw a - b 71 AT sw + b 72 A p R 
From (7) 



(1.29) 



A w = — A P R _ £i± A P ££1 A 



& a 23 w Cl 23 ^ a^3 -/ & 



% 



(1-30 ) 
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Substituting A P s from (2 3) and Aj* z from (20) into (1.30) 




or, 

AN r - AP ft - b 74 AT s - b 73 Aj> s ( 1 - 31 ) 



From (18) 



AN. . £ii 4X V + AP S + — AT S (I ' 32) 

s " a 33 V a ^3 a 33 b 



Substituting A P s from (23) into (1.32) 



AW, = A X v + f a »j a »l + fiiW + AP_ 

a 93 [ a 3 3 a 9 3 J a 9 3 ■' 



or, 

Aw s = b 7e AX V + b 77 At s f b 73 Ayo s ( i . 3 s ) 

.Substituting AW & from (1.31) and A W^ from (1.33) into 
(1.1) and rearranging ; 

4A = ( b 2- b 7«" ^r^s) A/ s + ( b 2 b 77 - b £ .t 74 ) AT S 

+ ( b f b 73 ) A P B + C b 2 b 76 ) AX V 

or, 

V« * C ‘i 4 A + C 12 AT i * A p B + D „ A *v (a) 

Substituting AW ft from (1.31) into (1.4) and rearranging ; 
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«© 

AT S * b 4 • b 7 5 A A + ( b 3 + b 4 . b 74 ) AT 4 - b 3 AT Sw - b 4 b 73 A P B 



Qr, 



A = ^ 21 A p, g •+ ^22 ATg + C23 ATg w + C 23 



A Po, 



C B) 



Substituting AW & from (1.31) into (1.11) and rearranging 

AT sw — - b is , b 75 A - ( b 4 g . b 74 + b j_(3 ) AT^ ■+ b iB AT S 



1 sw 



+ bis . b 73 A P ft + b^ 2 A W| + b i4 A W a 



or. 



0 

^"^"sw 3 ^ 31 A^g + AT s + C 33 AT sw + C 3 gAP B 

+ + d 33 4w a (c!) 

Substituting AW R from (1.31) into (1.20) and rearranging ’ 

A * = b 4 2 : b 7a A /s + b 42' b 7 4 AT S + b 38 A X - b fl7 A W 
+ b 39 Aw vv - b|3 AT bw + ( b 36 - b 4.2 • b 73 ) A P & 

* b fo At w + b +i A W [ 



or. 



AX 



« ^4-1 A /s +^42 AT s + d 44AX + CJ 45 AW + d 46 AW w 



+ C 47 AT BW +^48 AP 6 + C 49 AT W + *44 AW i 



(D) 
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Substituting AW 8 from (1.31) into (1.25) and rearranging ; 
AW » - ^66- Ws ^/s “ • ^74 At s - b 63 AX + b 62 AW+ b 60 AW w 

+ ^67 at bw + ( b 66 ‘ b 73 + b 61 ) A P A - b 6+ AT w + b 6 ~ AW{ 

or, 

o 

AW = e 8l Ayo + c 52 AT 6 + c; s4 AX + £ 55 AW + C 56 AW w 

*£57 AT bvv + £ sa A P ft + d 59 A 7 W + C>s4 AW- ( £ ) 

Substituting A W R from (1.31) into (1.22) and rearranging • 

* 

AW w 3 - b 50 . b 73 A/? s _ b 5j . b 74 AT S - b 46 AX + b 4s A W 
+ b 47 A W w + b sl AT &W + ( b so . b 73 - b 4+ ) 4 

- AT W + . b 43 A 

or, 

* 

^ Ww " '^ 6 i A / S 4 d 62 AT S + ^64 A >< + ? 6 S AW + C 66 AW w 

+ d 6 7 AT RW + d68 A P R + C e g AT W + D 64 AW' (p) 

From (1.29) 

AT bw . . b 7j AT S „ t AT bw + b 72 A P„ + b fc9 4 *, + b 70 AW a 
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or, 

at 6 „ . c 7a AT S „ + C 77 AT m + e 78 AP 6 + D ? 2 AW }+ D„AW a (<J) 

Substituting AW & from (1.31) into (1.17) and rearranging • 

^ “ ^33' ^73 A y°s, - ^35 ' ^74 AT S + b^ AX + b ft0 Aw 

' b 32 AW w + ( b za + b 3g . b 73 ) AP fi+ b 3a AT w -b a+ Aw. 

or, 

© 

APr * + ^82 AT s + d 84 AX + (1 <J 5 Aw + c %& Aw w 

+ c '88 ^ + ^89 AT W + ^84 AW i ( H ) 

Prom (1. 15) 

0 

AT W = b i5 A X + b £ ^ A W + b^ 6 A W w + b^ A + b, 7 AT W 4 ^ AW; 

or, 

AT W m (-<)4 A X + dg5 AW + C 3il A W w + dgg AP^ + d <3 g AT W + D 54 AW - ( X ) 

From (1.13) 

1 

Ay = b l9 AX + AW + b i0 AW w - b i7 AP B + b 21 AT W + b 2 . 2 AW( 
or. 

Ay . ^ 4 AX * i ..= AW + c: m iw » tC; «* A V <; .. 3 4 V l i,f w i ( J ) 
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The state variable equations from (A; to (J) were rewritten 
in matrix notation as shown in Figure 4. The state equations 
were solved for various percent step changes to input 
variables ( AX y , A W p , Aw j ) / using the IBM simulation language 
CSMP - IH as shown in the next section. 
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> 

X 


u_ 


< 

> 


3T 


<1 


<1 


<J 


<1 



1 

o 


o 


o 




o 10 


o t0 


o 




■* 

d* 


~xn 

o 

or 4 


o 


o 


fO 

o ' 0 


o 


o 


o 


f0 

cT 


o 


o 


o 


O 


o 


d” 


o 


o 


o 


<N 

cT 


o 


o 


o 


1 . 


o 


o 


o 


o 


o 


o 


o 


o 


o 

1 


*4" 


1 

i- -•?■■■ 


<1 


d 5 


X 

< 


> 

<1 


3? 

<1 


3: 

<J 


of 

<3 


h- ? 

<1 


^ 1 
< . 




1 

o 


o 


o 


o 


o 


o 


o 


o 


o 


1 

o 


o 


o 


o 


CD 

•* 

TJ 


cn 

10 

U 


a> 

vO 

TJ 


o 


o? 

cO 

XJ 


oy 

ay 

Xj 


cry 

o 

Tjf 


cO 

■*4 


00 

C4 

O 


00 

c9 

TJ 


CO 

XJ 


CO 

10 

■o 


00 

vO 

U 


CO 

N 

XJ 


cO 

flO 

xj 


cO 

ay 

XJ 


cO 

© 

TJ 


o 


o 


O 


h* 

XJ 


10 

XJ 


vX> 

■o 


N 

N 

XJ 


O 


o 


O 


o 


o 


o 


sO 

XJ 


U5 

10 

XJ 


so 

Vp 

o 


o 




40 

■u* 


*o 

o 


o 


o 


o 


10 


10 

10 

-o 


ID 

40 

XJ 


o 


to 

cO 

XJ 


10 

Xj* 


10 
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Figure 4 - STATE VARIABLE MATRIX 



III. COMPUTER RESULTS 



From Appendix A , it is seen that the a coefficients are 
in terms of steady state values. These values were 
determined from the data in Whalleys thesis (2) which are 
repeated in Appendix E. A Fortran IV program to find the 
a,b,C,D coefficients which appear in the CSMP program is 
shown in Appendix F 3 and the calculated values are shown in 
Appendix G. Appendix H is the IBM simulation language CSMP 
program that was developed. Only input variables were 
changed for each computer run. For example, results for a 554 
step change of A X v , a 10% step change of AVljand a 10% step 
change of Aw t were calculated. The output response curves 
were plotted using the NPS CALCOMP PLOTTER. 
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A. TRANSIENT RESPONSES FOLLOWING A STEP CHANGE IN THROTTLE 
SETTING OF 5* 



At a particular steady state operating condition the 
throttle valve is suddenly opened while still keeping the 
previous air and fuel flow rates and feedwater flow rate 
constant. Because of bubble formation in the steam drum 
liquid, the steam drum water level (FIGURE 9) swells for 
about 40 seconds and then shows a steady falling in level at 
higher flow rates of steam. The steam drum pressure (FIGURE 
3) immediately begins to fall. Following these sudden 
increases, the steam mass-flow rates from the steam drum and 
the superheater (FIGURES 10 and 11) show slight declines 
because of the effect of the decrease in steam drum 
pressure. The superheater outlet pressure (FIGURE 10) acts 
in a similar fashion to the steam drum pressure, that is ,it 
shows a steady decline to a new steady state condition but 
it is slightly lower than the steam drum pressure because of 
system pressure drops. The superheated steam temperature and 
superheater wall temperature (FIGURES 5 and 6) both decline, 
although the latter temperature is slightly lower. The 
riser tube wall temperature (FIGURE 8) and steam drum liquid 
temperature (FIGURE 9) also decrease due tc the throttle 
change. The superheated steam density (FIGURE 5) shows a 
rising characteristic to counteract the loss in superheater 
pressure and temperature and the gain in steam mass-flow 
rate from the steam drum (FIGURE 11). The riser and 
downcomer flow rates (FIGURE 7) are reduced to a lower value 
but steam quality (FIGURE 6) is increased to a higher value 
than previously to preserve the energy balance. 
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Figure 5 - SUPERHEATED STEAM DENSITY C AyO $ ) V S TIME AND 
SUPE8HEATED STEAM TEMPERATURE < AT S > VS TIME 
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Figure 9 - STEAM DRUM LIQUID TEMPERATURE < AT W ) VS TIME 
AND STEAM DRUM WATER LEVEL < Ay ) vs TIME 



57 



G.cn 



tz.no 



[Xin J 1 
te.na 



24. an 



30.Q!ki 




txj a j 1 

-O.OD S.C1D 12. CD 18- DC 24-Du 30. 




-0.00 B.O 0 


12. QQ 


ift.na 


24. on 




30. BO 


TIflE 




CX1D ] 1 








Figure 10 - SUPERHEATER OUTLET 


PRESSURE 


< AP $ ) 


vs 


TIME 


AND STEAM FLOW RATE AT 


SUPERHEATER OUTLET 


C Aw s ) 


vs 


TIME 






58 



1-92 i-3a 2-04- 2.JQ° “-9.DD -6* Dll -3.DD -0.0 

DW5 KID J ~ 2 DPS WO 1 1 



CO'H CQ-9T 00-9 CO '0- 

_ [ cruj 9 kq 



[X1D 1 




Figure 11 - STEAM MASS-FLOW HATE FHOM STEAM DRUM TO 
SUPERHEATER < > V S TIME 
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B. TRANSIENT RESPONSES FOLLOWING A STEP CHANGE OF 10% IN 
FUEL FLOW RATE. 



The effects due to a change in the fuel flow rate are 
less than those due to a change in the throttle valve 
setting. However, the system requires a longer period of time 
to reach a new steady state condition. Because the 
combustion rate is increased, the evaporation rate is also 
increased making both the steam mass-flow rate at the 
superheater (FIGURE 17) and the steam mass-flow rate from 
the steam drum (FIGURE 18) increase in a similar manner. 
Initially, density in the risers falls quickly in response 
to the increased firing rate which causes a rise in 
superheated steam density (FIGURE 12) . An increase in steam 
flow from the steam drum without an increase in the 
feedwater flow rate results in a drop in steam drum water 
level (FIGURE 16). The steam arum pressure (FIGURE 15) and 
superheater outlet pressure (FIGURE 17) rise monotonically 
as time increases. Also superheated steam temperature 

(FIGURE 12) , superheater wall temperature (FIGURE 13) , 

riser tube wall temperature (FIGURE 15) and steam drum 
liquid temperature (FIGURE 16) all increase with increasing 
fuel flow rate. The abrupt fall in riser mixture flow rate 
(FIGURE 14) is probably due to the sudden drop in the level 
of the initial density change in the riser. 
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Figure 13 - SUPERHEATER WALL TEMPERATURE l AT SW ) vs TIME 
AND QUALITY OF MIXTURE LEAVING RISER C AX ) VS TIME 
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Figure 16 - STEAM DRUM LIQUID TEMPERATURE C AT w > vs TIME 
AND STEAM DRUM HATER LEVEL < Ay > VS TIME 
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C. TRANSIENT RESPONSES FOLLOWING A STEP CHANGE OF 10% IN 
FEEDWATER FLOW-RATE 



The main effect of increasing feedwater flow rate while 
maintaining throttle setting and air-fuel flow rate constant 
is an increase in the steam drum water level (FIGURE23). 
Because feedwater is added into the steam drum, steam drum 
pressure (FIGURE22) suffers a small overshoot followed by a 
decline. As a result, a drop in both steam mass-flow rate 
from the steam drum (FIGURE25) and steam mass-flow rate at 
the superheater outlet (FIGURE24) occurs. Flow around the 
riser-downcomer loop (FI3URE21) is only slightly affected. 
Since the firing rate is unalter ed, the riser tube wall 
temperature {FIGURE 22) and the steam drum liquid 
temperature (FIGURE23) decline. The superheated steam 
temperature (FIGURE19) and the superheater wall temperature 
(FIGURE20) rise because of a drop in steam mass-flow rate 
from the steam drum. The quality of steam (FIGURE20) shows a 
general lowering. The superheated steam density (FIGURE19) 
shows a decline which is similar in form to the drop in the 
superheater outlet pressure. 
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Figure 19 - SUPERHEATED STEAM DENSITY < Ay° s > VS TIME AND 
SUPERHEATED STEAM TEMPERATURE < AT S > VS TIME 
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Figure 20 - SUPERHEATER WALL TEMPERATURE < AT $W ) vs TIME 
AND QUALITY OF MIXTURE LEAVING RISER C AX > VS TIME 
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Figure 21 - RISER MIXTURE FLOW RATE C AW > VS TIME AND 
DOWNCOMER LIQUID FLOW RATE < AW w ) VS TIME 
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Figure 22 - RISER TUBE WALL TEMPERATURE < AT QW > VS TIME 
AND STEAM DRUM PRESSURE < AP & ) VS TIME 
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Figure 23 - STEAM DRUM LIQUID TEMPERATURE < AT W > VS TIME 
AND STEAM DRUM WATER LEVEL < Ay > VS TIME 
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Figure 25 - STEAM MASS-FLOW RATE FROM STEAM DRUM TO 
SUPERHEATER < AW B > VS TIME 



75 



•0.00 -4.00 4 00 

DVB [X10 1 " 4 



IV 



CONCLUSIONS 



This analysis is only the initial step toward a better 
understanding of a naval steam generating plant both from 
the dynamic and from the control points of view. The 
analysis concerns a boiler system which is divided into four 
parts. The major difficulty in the analysis is the fact 
that the whole system is very complex and contains numerous 
variables which are extremely unwieldy to manipulate. Large 
numbers of relationships are non-linear so it is necessary 
to generate a linear form by the methods of perturbation 
theory. After setting up the state variable equations, the 
IBM simulation language CSMP-HT was used to solve for the 
open loop transient response characteristics of the state 
variables due to variations in input variables. Throughout 
these calcu lations, all the variables appearing in the 
transient response curves are not absolute values. 
Hather,they represent a small incremental change from some 
particular steady state operating condition. Also, it is to 
be noted that trends in some of these results differ 
markedly in form from those appearing in Whalley (2) . Future 
studies will be undertaken to determine the causes of these 
differences. Meanwhile, t he benefit one may derive from this 
basic analysis is its use in the design and construction of 
a new type of boiler controller, that is, a multivariable 
control system. 
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A 34- 


\ a 35 


a 35 


( ^■3 2*^55 i 


a a6 


V a 35 


a 3 5 


a .^- b S6 




a 35 




a 3 S.' b 37 




35 




a .12- b 38 




a 35 




a 32*^5 9 




a 35 





88 



6a 



6 9 



•( 



a 08 Q 5o 

fl 03’ a St 
Q 40 ’ ^-4a 

^38- a 62 



+ 



71 



72 



73 



74 



75 



\ <140 • ^4.5 

_ Q 38- a 49< b 7 

a 40’ a 4fc 

_ ^33- a 53 ' ^37 

a 4c • Q-ss 

_ / Q 26 _ a as a 3a ^ 
V, O 23 q . 23 / 

°-ZA- q -16o 



a 39 ' a 54 
a 40' °-53 



a a8 a 49- Pfe 
a 40 ‘ < 3 - 4-3 

a 3 9 a 48‘ ^8 

a 4o-a 49 



tt 24 Q -101 



76 



7 7 



u 94 

0-9J 



■( 



q 95 -Q-ioq - 0.-36 

a 93 093 



b 

7 8 



Q 95-Q-101 

a 93 
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APPENDIX C 



C 111 TR IX COEFFICIENTS 



c 




(b . b -b . b 


1 1 




2 78 1 75 


c 


= 


(b .b -b .b 


1 2 




2 77 1 74 


c 


— 


b .b 


i a 




1 73 


C 




b .b 


2 1 




4 75 


c 


— 


(b +b .b ) 


22 




' 3 4 74 


c 


— 


-b 


2 3 




5 


c 


- 


-b .b 


2 a 




4 7 3 


C 


— 


-t ’ .b 


3 1 




IS 75 


c 




- (b . b +b ) 


32 




IS 74 16 


c 


— 


b 


33 




13 


c 


= 


b .b 


36 




IS 73 


c 


- 


b .b 


4 1 




42 75 


c 


— 


b .b 


42 




42 74 


c 


= 


b 


4 4 




38 


c 


= 


-b 


4 5 




3 7 


c 




b 


46 




39 


c 


— 


-b 


4 7 




43 


c 


= 


(b - b . b ) 


48 




36 42 73 
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c 


= 


b 


49 




40 


c 


— 


-b .b 


SI 




66 7 S 


c 


- 


-b .b 


52 




66 7 4 


c 


= 


-b 


54 




63 


c 


— 


b 


55 




62 


c 


= 


b 


56 




60 


c 


— 


b 


5 7 




67 


c 


— 


(b . b +b 


58 




' 66 73 


c 


= 


- b 


5 9 




64 


c 


— 


-b .b 


61 




SO 7 S 


c 


— 


- b . b 


62 




SO 74 


c 


— 


-b 


64 




46 


c 


— 


b 


65 




4 S 


c 


— 


b 


66 




47 


c 


— 


b 


67 




S 1 


c 


- 


(b . b -l 


68 




SO 7 3 


c 




s -b 


6 9 




48 


c 




: -b 


73 




7 1 



c = b 

7 7 68 

C = b 

C 

c 
c 



7 0 


72 






= -t . 


b 


8 1 


35 


7 S 




II 

1 

rr 

• 


b 


02 


35 


7 4 




= b 




04 


3 1 





) 
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c 

85 



c 

86 



c 

88 



C 

89 



c 

94 



c 

95 



c 

96 



c 

98 



c 

99 



c 

10 4 



c 

10 5 

c 

10 6 



c 

10 8 



c 

10 9 



30 

-b 

32 

(b +b • b ) 

29 3 5 73 



33 



25 



24 



26 



23 



27 



19 



1 8 



20 

-b 



1 7 



21 
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APPENDIX D 



D MATRIX COEFFIENTS 



D 

1 1 



D 

32 

D 

33 



D 

44 



D 

54 

D 

64 



D 

72 



D 

73 



D 

84 



D 

9 4 



B 

10 4 



b . b 

2 76 

b 

12 

b 

1 4 

b 

41 

b 

65 

b 

49 

b 

69 

b 

70 

-b 

34 

b 

2 3 

b 

22 
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APPENDIX E 



LIST OF CONSTANTS 



L 

S 




45.0 


A 

S 


= 


0.239 


A 

c 

ps 


: 


0.5 

0.53 


h 

S 


= 


1390.0 


h 

a 


= 


1204.9 


w 

s 


= 


4.0 


M 

S 


= 


1900.0 


C 

S 


= 


0.23 


K 

gs 


— 


0.94633 


W 

f 


— 


0.32 


T 

gs 


= 


2800.0 


T 

sw 


= 


1250.0 


K 

as 


= 


0.0 


w 

ac 


= 


7.0 


w 

a 


= 


5.3 


K 

s 


= 


4.8848 


T 


— 


1200.0 



S 
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4.0 



W 

B 

C =0.37 
C 

T = 3200.0 
C 

f = 68.22 
S 

jO & = 0.758 

A = 2.95 
r 

L = 10.0 
r 

p = 16.52 

A = 52.6 

X = 0.031933 

K = 0.00001 
B 

C = 0.007637 

W = 178.97 

f = 0.76095 

r 

D =0. 146 
r 

g = 32.2 

W = 178.97 

H 

M = 12000.0 

B 

C =0.28 
B 

K = 0.00406 
C 

h = 390.0 

fg 

h =0.0 

UD 

h = 4 30.39 

K = 3.6632 

gB 

K = 7.5E-10 
r 

T = 3000.0 
gB 



95 



930.0 



T 

bw 

K =0.0 
ar 

K = 0.059806 

g 

T = 890.0 
B 

f = 2.7014 
D 

L = 10.0 

D 

D =0.3 
D 

A = 0.426 
D 

V =40.0 
B 

A = 26.4 

K = -0.39157 

e 

T = 740.0 

i 

T = 894.38 

H 

H = -1.7742 

e 

h = 405.0 

BW 

M = 1680.0 

K = 0.003 

T 

K = 500.0 

f 
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APPENDIX F 



FORTRAN IV PROGRAM 



- r-f« t tor*** — (-/v^v - ) 

l S = ^ « 0 

AS = 0,2 3 ■; 

RHO S = C • 5 
CPS = 0 .53 

hS=l3frO . C 

KB =1204 • 5 
W S = 4 • 0 

-fr S='LVOO vO — — — — 

C S = 0 • 2 3 
h r; s = 0 . 4 4 5 b 3 

TTF^TTI/ 

TGS=2300 .0 

■XS.W-=-l-2-5i^ 

K A S = 3 • 0 
WAC= / . 0 

WA = 5; - 2— — 

.< S = *+ • 6 •' 4 : 

TS = I 1.2 00 .0 

*1 0 = 4# O 

C C = 0 • 3 7 

■Jt: =33racu-c 

F S = 6 3 . 2 3 
PH0ri = 0 . 7 5 m 

S-j 

LR = l ) . 0 

rHO lo . b 2 

k h f J .n - ;> ? • t 
X=0. J31933 

Jca-G.occci . — - - 

C=0. J07627 
W = 17 2 . 37 

F k= 0. 70033 — — 

0 k = 0 . 14 5 

22.2 

777^TT77T7 

M 5=12 3-J0 . J 

lOa=2U2-B 

ftC =u . uO 't Go 
HP G =23 0*0 

-h7/~nP=Gr^9 

h = 4 3 0 . 3 J 
K 03=3 .6632 

T^=77“'J7-n 

TGb=3 J0C. J 

=0 3 ,0^0 

k*\R= ). 0 
KG=0. 0 5<= c06 

Ii^3 ;n.Q 

F J = 7 . 701 4 

1 L) = l ) • • J 

UTT=7 .7 

A 0 = 0 . 4 .7 6 

_V_d_ g_4_.3_«.Q 

A = 2p • -f 
K c =- J. SS l? 7 

r /i = ? ^4 .3 b _ 

t*i t = — l • ' ( i 2 
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to =1630. U 
K T = J . 302 



=5 00 . C 
A l = l . 0 
- A 2 = - l »Q 



A : = L S A S 

A4 = l S*A S^ishOS* CPS 
A 5 =K3 — n o 

A 6 = - 1. . 0 

. A7a-.TiSttr.PS 



A3 =1.0 
A c , = — l • 0 
■ 10 = MS*CS 



a n = i. . o 

A 12= (6. 6* < GS ) / ( 
AI3=K'GS*vsF**Q .6 



■ Ftttt.j . A) *( TGS-tS'a) 



-st^t 



'.AS"/v;A tP 






. T - V t • C 

A 1 5=1.0 

A 1 6 — ( 0 . 'i ^ C S / v': r' ^ 0 ♦ 2 ) f T S Vi ~ T 3 i 



i 1 = K S — a 3 tt ” J . c 
AlS = UO/ (2.0*CC) 
4, 



A 20= T C- T GS 
A 2 1 =A2 0 

~A~ 2 2 = •“ ( IvF f vi A i 

A23=2 .0*FS*'/.3/RHC8 
A?4 = l 

A 2 5 = — ( r ) 

A2e=i .o 

?* ( l . 0 ARH C V. - 1 ■ . 0 /•* hO P.J- 



A2 3= ) .5 A A < *!. KttXtt.K b * ( ««-r?tttt 2/ RHG P**2 ) 
A29=0.5ttCtt(RhG + RI-Cv> ) 

A A 0 — 1 • c 
A3l=- l .0 
A 3 2= L . C 



A^=-L.O 

A3 4 = 2.C*«*(FRttLR/nF; + L.5 )/( G*A ? **2*RhD ) 
_A?4L = Li U-i =iA^AuU 



A3 6= ( 1.3 /RK13- l . C / 1< r C V ) *( ( W** 2 / ( f-tt A Rtt* 2 ) ) * <Fft*LR/CR+L.51- 
1 ( j?H0tt«2ttL^) ) 

- - A 3 - 7= - .rvG - ->y. rrHrr A A ft * * 2 * F C r* ) 

A 38= l .0 
A3 9=- l .0 

A ■+ J= v ! ti tt l y 



A'+l = D.5*ARttlj^KCtt ( RFJ + PK1 4 ) 

- A A 2 = 0 . 5 v A K -r i_ .9 ai h F G tt ( 5442 iS. 5.0 4 ) 



A 4 2 = - 1 . 2 
A +4 = !-iv< C- h 

A'+ 5 = -' , ■ * KC 

AA6=-/ittl-Po 
A + 7 = ^ 

A 43 =1.0 

A49 = ( -KGOtt^r- .ft ) - ( 4.0ttKR*TG3** ? ) 

in = ( -K „ F4-V!) . A ) - ( 't -OttK A 4 T .-> * *?% 2 ) 



A 5 l = 0 . 6 *K o 3* ( T G a - 1 8 * ) / * F*tt ) . A 
A3 2=2.C*KAR-s (V. AC-ir.A ) /VAC** 3 
- y r5- 2 = l -i- e - 



A 5 4 = 3 . 3 tt K 7 * ( T o V\ - T ls ) •* * 2 • 0 
A 5 5 = - A 5 + ' 

A5c-i . O 
A57=- l .0 

_t\5 3 =2 - P* r'i-V G/l-CfO .5) / ( AO**? ) 



A5C = i_l/ ( GttA'i ) 
AoO-Vli^K H. 

~ A 6 t — — P Ft _ -i * A 



Ac2 = K'- 
Ac 5 “/ 



AC ■+ = A 

Ao 5=- JL • ) 
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flfeft=TY.»t 1 . j-X) 
Aft7=W* ( .1 . J-X ) 
A6S=-W*7B 

A6 »^— T- 

A 7 0 = — vn 
A7J.=-TW 

A72=-WE*s<C 

A?2=-HBIn*Kc 

A 7 4 = A 

A75=A*RHCw*T a 
A 7 6=— K c 

-A 7^ = ^4AM^ 

A 8 S = 1 . G-X 
Ac 

A 3U=-l . u 
A 9 l = - < c 

A9 2 = 1 . 0 

A93 = l . 0 
A 9 4= 0 . 4 2 
4 W - S == 0. -3 099 - 1 - 9 — 
A 9 6 = — 0 . 0000 )5 
A 9 7 = K T 



^ ^ u 

AO v-XC: 

A 100=24 0.0 

A l J l = 9 5 9 C J . J 

A 1 0 ? = < F 
B1.=A1/A3 

r. 2=A 2 A A- 3 

B3= ( A 6* A 1 7 ) / ( A 4* A 1 9 ) + ( A7/ ,i4 ) 
b4 = ( A 6* A lo ) / ( A4-A.15 ) M A 5 / A 4 ) 
d j - { A t * A 1. 7 i / ( <-. 4 ^ A i 'j ) 

uc=( ( A 19* A 1.0 2/ A 1.0 )- ( A1.2/AU ) + ( 7 70 /A 1.5 ) ) /< ( A 1.3/ A l U-( A 27 /A 1.0 ) ) 
8 ? = ( A 1 3 / A l l ) / ( ( A l 2 / A U ) - ( A ? 2 / A 1. 5 ) ) 

be = ( A l 4- ( A rr / A 1 0 ) ) / ( l A 13 7TT I. ) -TJTZ7TT3T1 

b 5 = ( A 1 2 / A L l ) + ( A 1 3 * H 8 / A 1 1 ) 

c.LO= U 7 » A U / A t i,)-.( A l 0/91 U 

BL I. =A 1. 4 - { U.*A 12 /A L U 
rtl2=A3*09/Al0 

8 1 3 = ( A 0 " i 1 . 0 / A 10H ( ( A i v T r 70 - / - ( - A+e^i -5~ H 

Bl4 = A3*i> ll/A 19 

815= ( A^ALc) / ( A 1 0 * A 15 ) 

b 1.3 = ( A 9 * 'A 1 /l/IOUTTm 

b 1 7= A 7 6 * A 9 7 / A 7 7 

UVS=ARH/A7.7 

bl5 = Al>*7 / A77 
B 2 0 = A "! 0 / A 7 7 

ti2 1 -= Av 1/ A? 7 

t>2 2=A C 2 /A 7 7 

323 = {A75*317/m74)+(A4'»*A6 7/A74)-(A?7*A7 3/A74)MA72/A74) 

ti'-' 4= Ue c/4 7 t)-( 1-VA74J 

b25=( A 6 8 /A 74 J- ( A 7 5 L9/A74 ) 

ZJLL= LA-7 1 / A 7 4 3 - ( O--7_5702.iq/J\72tJ 

b2 7 = ( A 7 0 + A 73 )/A7 •>- ( A 7 5 * 8 7 l / A 7 4 ) 

623= ( A65/*74) - ( A? 5-822/A -4 ) 

r~29= ( Afc1.t.jl7/4CC) — B Ao2* VJ7/A 6 0 ) 

B30=( A63/A6J) - ( A 6 l * R l 5 / A 5 0 ) 
b2 1 = ( Ao4/«6 0)- ( A 6 L*B L r /A60 ) 
bb7 = Ao i. - r ' 0/ AoO 
•3 3 5 = ( A 6 2 / A 6 0 ) - ( Ao i* P. 2 l / Aft )) 

f . ^.4 =Aba_* b 2 2 XAaJ3 

B55=Ab5/AO J 

r 3 6 = ( A ^ 4 / A 4 ? ) - ( A 4 l * 1 1 2 r - / A 4 2 ) - ( A 4 3 * A 5 5 * A 9 ^ M 4 2 / A 5 2 ) 

u 3 7 = A 4 l l -< u 0 C /-~V4-g 

r 3 8= ( A 4 5 / A- 7 ) - ( A4 l * B U / A4 ) 
r 2 9 = ( *w 4 /A 42 ) B ( A -* 1-:! 2 "'/A 4'’ i 

l,4.: = i Ah : . / -\4 1 I * 

6*+ 1 = h 4 l / A -t 2 •’ ;> 2 4 

642 =A4 1 Bi -i9/A4~ 
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843 = ( A43*A54)/ l A 't 2 * A 5 3 ) 

644=( A2 /*636/A29) + ( A? 2*B29/ A29 ) 

- B 4 5 » < *•£ 7 »->? 7-/* - gO ) - t A 2 6 ♦ 35 0 -/ 4V-?9-) -* 4 A 3%/ A 29 ) 



64 6 
647 
o4b 
649 = 
B 6 0 = 



6= ( A 2 7 * •-> 3 o / A 2 9 ) M A2 **‘33 1 / A 7 9 ) 

7-( A 2 3*322/A29)+( A5C/A29)-( A?7*B2S/< 

i ™ _ J ^ . p j i p < * — ■ -rrr ^ * ■■■ — r a f iu n * ^ 1 



TKTTT 



T7V27”6 40/42"' ) + ( A"Z F "* B ~3 
( A2«*334/A29 ) - ( A2 7*3*1 /A 29 ) 
( a 7 7*34*/ A 29 i-( L2Z*h25JA2 2 1 



/A 29 ) 



tj5l = A2 7*8^3/ A29 
B 5 2= ( A 5 3 / A 5 7 ) + ( A 5 9 *U4 7 / A5~ ) 
653= t Ad r - / A 3 7 j * ( A 5 9 * 6 i A / A 5 7 ) 



= A 5 3 * 6 -t 5/ A3 7 
= m 5 9 * 6 4 6 / A 5 7 



B i> 4 = 
h 5 5 = 

3 3 c = A o •> c + o / A 5 7 
B57=A69*o49/A57 
65B=45?*03 Q/a37 



: A 59*6 5 1/ A 5 / 

= ( A3 2*6 5 2/ A ?. 5 ) - ( A? 7 / 42 5 ) 



66 9 
B 6 3 

-f*Hr . . . 

66 2= ( A3 2*6 5 4/ A 3 5 ) - ( A 3 4 / A3 3 ) 
P 6 3= ( A 3 2 * 3 6 5 / A 3 5 i + ( A 3 6 / A 2 5 ) 



= ( - A3 ^ 7-V - : - >) — t-A-3- 



;-§4M- 



4- 



33o/ 43 5 



no- 

= A 3 ?. * B 5 7 / A 7 '= 



bo *t= 

B6 5 = 

B6c=A32* r i5 6/A2 3 



6 6 7 

d 6 6 
-&-6 



=A32*c59/ A .; 



= ( A 3 
=+- 



A 5 0 ) / ( A 4 0* A 4 7 ) + ( 

■ ^ -.‘f 



: A 5 19 / ( - V 






)-( A - 3 



3 C *A54) / ( A40*A £ 

\ /. c i 



) 



^ 6-) - / < ^ -4 G*A 4- 3 )- 



D 7 0 = ( A 3 3 * a 5 3 I / ( A40*A4- ) * ( A? S* A 4 9*88 ) / ( A40*A49 ) 
671= " ' ' 



4Ui 

n . ^ 
4- 
b7S> = 
“TTTt" 



( A 

( A 



=t-*T 






* A 5 1 r 






“ b 7 ) / ( A 4 0 * A 4 1 ) 
s 7 1 / ( M0*A5^ 






4 2 4* A LU0/A2 3 
A 24* A l 01 /A 2 3 
A ^ 4 1 A 9 2 



e 7 7= ( A 9 3 *A 1 . 00 / A93 ) + ( A^6/A93 ) 
6 7S=A9 5* A l.61/A . 9 . 'j .. 



LI 1= ( 62*676) -( 6 1*37 5 ) 

Ll2 = ( 6 2 * H 7 7 ) - ( 01*674) 



C? 1=34*6 75 

L 7 6 = 3 3 + ( -4*0 / 4 ) 



TTT^ 
L2 6 : 
-C4L= 



=-b4*B7: 

1 ~ 6 1 6 * n 



13/.= 
C 3 

trf 



- ( 3 l j ■= 6 7 4 ) - 6 l 6 
6 13 

3 L 2*;> f 7 



C4l = 
C 42 : 



642*075 
■342*6 / 4 



^T3T 
= -6 3 7 
:| ) 3 3 .. 



L4 4= 
C45 
X-46 



C47 = 
(.46 = 
-tr4*4 



- o -t 3 

6bo— (6 42*67 5 ) 
tHf 4 



Co) =-665*675 
C5? = 

T54^ 



=- oo ■: 



*6 74 



C 5 5 
-C2L6 



=-bc 
= 3 o 2 
=3-2J- 



C 5 7 = 
C5< 



-644- 
C 6 l = 
C o2 = 



oc 7 

( .666*673 ) +Re l 

■-t rc - 4 

; - B 5 0 * 6 7 5 

-650*3 74 

- '4 
=o45 

ihai 



of: <t= 

Co 5 = 



C6 7 = 
C6o = 
0 6 °~= 
r.73 = 



c5 1 

( 650*3 73 ) — 6 4 4 



TT 

-67 l 



100 



C 7 7=06 9 
C / l=o TT 



C 3 1 =-l 35*6 



, 5 * 02 - 



C3 4=8 i L 
C 8 5 = 13 3 0 
C. .10 - 



C38=B29+ ( o 3 5 * b 73 ) 
C o 5=B 3 3 



L‘> "+=0 Z o 
C 95 = 3 ?^ 



C v 8= 62 i 
099=82 7 

■ - a i 






C l 0 3 = 8 1 3 
C 10 fc>=P 20 



l. i. 0 3 =— o ]. 7 
ClU 9 = 3 2 L 



f l l=o 



t: 



2 2 =b 1 2 

133=3 L 4 
■ L. 4 4 = 3 4 1 - 



034=365 
0 64=34 5 



"777~2"="3T79"“ 

073=370 




5 TT = , t r 3 Tr 75 TT 5 TT C 8 = =" . 1 : 1.3 
5 >- C 6 2 =,tl5.5,5X, 5HC64 =,515.6/) 

r r jj ,r f 7 .r. 43 ,r 69 t r. 7 ~.C 7 ? , C . ZJ 1 m!£ l-,JCJi 2 -,C- 34 ., L.- 25 . , 

1 0 86 , C 3 1’ , Co 9 , C6 4 , C 9 5 ,C 5 6 »C 9 3 , C9 9 , C 10 4 , C.$. 05 , Cl 06 , C 1 06 *C 1C 9 
200 FuRMAl ( 1 1- , 0 X , 5 H 0 6 5 = , F 1. 5 . 3 , 5X , 5HC6 6 = , E 1 5 . 8 , 5 X , 5hC6 7 .3// ♦ 

5 a , 9* ~ C o5 — = * C 1 5 . 4 ,-3 X , 5 F 0 7 _ : 



4 X , 33C o 

99X , 5 H C 6 1 
• t< .,T: 2 14 



A l \ - 



LQill. 



- 1 - 9 X-, 



r , 8 * 1 0 o ~ — =~r 
2 9 X , 5 rl C 7 7 =,8.16, 
59X , 5HC52 = ,FL5 



-V5r 



F L V . - 5 -/-/- 



4 *3 X , 4 = L o = , l 1 6 • " , 

5 5 X * 5 H C 5 4 =,315. 5 , 
.4 SLa , 3H(. 9 6 — = , F I. 6 » 3 > 



, 5X,5HC7= =, r l 5.8,5X,5HC31 =,E15.8//, 
3 A , 5FC 34 =,F15.8,5X,5HC83 = f FJ5.8// t 
F 15 . 377",' 



TTn =: 



6X,3 l «0H =,!.. 15.3 ,5X ,5HCh' 

5 X , 5FC 95 = , F 1 5 . 8 , 5 X , 5HC 96 =,F15.8//, 
8 V T r,i-r’-,-i = t F 1 4.° , 5X , 5HT 1 n4= f F 1 5.P.// 



40 0 



7 5 X , 34C l O 1 - , E 1 o . o , 3 X 
5 5 X , 5-10 10 2= ,5 15 . c/ / ) 
— rrr -H"- 1 - — j.j) 9 ) - .l - " 



iFC l C 5 = , E 1 



5 X , 5 F- C 1 0 8 = , E15.R// 



F 0 KM A T ( I. F J ^.<,5 A HI l 



=,F15.o,5X 



tt 

680 



\. c - X , 5 H J 4 > : 1 1' 1 . • ‘ ,3- , 3 - 1 / '4 — , F 1 5 . 



J. A f 9 f t — y 

~7 t 77~TT^TT2 =7 

3 C . X , 5rlD c/ -t =i 



3 c Xt 5r 1 0 < 

jsin?. 



= /r l 5 • 



T~T5TT3 t ‘TT 
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APPENDIX G 
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APPENDIX H 



CSMP-HI PEOGRAM 



i£. 6 CrjNT: i MUOt ' S - SYSTEM '■COCLu'.G Pn,CGt» V i ^3 T-.h\-\.Si. AT Ol — OUT- b T 44 ? 

IN IT" AL 

CONSTANT T A l J I . = T . r TT^T ^ 

CONSTANT 871=0. 1.4 l 58- 32 , * 7 A = C . 2 ? r ‘ r + 0 1 , R r 5 =0 . 1 3 1. 50 +0 ' 

CONSTANT R 7 i. = 0 . A 2 u 0 r + 0 0 , R 7 7 = 0 . A 4 4 5 -- 0 ? , R 7 3 = 0 . 1 2052+0 1. 

CONSf AN T C ll = -0 . rTTTTTu > ,C !2=-0 .3 I. 4 2'-- J I ,'*.1.5 = 0. l.'M6"-03 
CONSTANT C 2 1=0 • 1 7 148+04. C 22=-C. U 098 + 0 1 ,C22 = 0. 5 1. 7 3 r 4- 0 + 

CQ.vsxant— gxx— - =■;> . 1 .o ■> i - j 1 . . c 0 1 1 a + c " , c : 2 _ .3 , 1 4 A ■ . r + a r 

CONSTANT C : 3=- 0 . 2 486E- C 1 ,C3 c = -0. 4 * -Of F -C2 »C4l=-0. :4;4'= + 31 
CONST ANT C 42= -U . be 75 f 1 - C2 , C 4 4 = -0 . a 7 + £ ^ + c 1. , C 4 5 = - C . U c 4~- 0 ? 

CONSTANT CV; - - J. 7 J/u-A-O' 7 , . C 4 f = j - ' V ' , ' 0=0 .- - 7 • 1 • • -Ob 

CONSTANT C 45= J. lC^tCh- 3 .1 ,C 5 1=- 3.4"- ) 9= + 2 6 O 5? = - 2 . I. .?4T= HO 4 
CONSTANT C 54 = -0 .63 3 + r +06 . C 5 8 =-0 . * 5 7 5 8 + C ? . 0.5 6 = 0 . c 7 +0 2 

CONSTANT c5 7 = 0. c- L 4- + J2 ,C58=0. Or g 4 - 01.1 0 = j . 1 . Ac '7 -+ 3 3- 
CONST AN T Cc l.=-0 . 70 xi E +0 5 , C 6 2 =-0 . i 7 A J r + C * ,C 64 = -r . 9 1? 5C + C5 

C ONSTANT Cu5 = -G . 7) ; V 4 . £ + GX , £ . A , - - . 0 . - ■, : ~ ~ £ 1 - Q 3 . C . . 7 = 0 - .. I 4 5 2 ~ t . , ? 

CONSTANT C6 5= 3. 5 5 J. (O - 0 l » C 6 ? = C . ; 1. 0 7 - + 0 ~ ,0 7 5 = 0.25 4 3 - - 0 2 
CONSTANT C 7 7 = - 0 . 8 S 7 0 8 • - C 1 , C 7 8 = 0 . 7 5 " c - 0 3 , C 5 l = 0 . 5 2 = r c + 0 0 

C :j ‘ 3 S T A N T C 2 2 = J . 1 3 J Tr~rtT^ ~ Tt rT T ! r^Tt : * * *r ■ - 1 ~ + 0 4 / ■ 7 - 0 . 1 I. 

CONSTANT C3i>= - J .2&C3c + C2 ,C3? = -0.7 4 378 + C 2 ,C 5 9 =-r . Sr, ARC + 03 

CONST ANT C94=fl . 464 5 ~ ±0 0 , C. 4 5 = -0 . j 5 ? 4 r - -0 0 <Fo = - Q . 00-0 + 

CONSTANT C 90= 0.6 5530-03 ,C99 = - 0. 2 ? 34 F- + 0 2 ,0.1'.* ■■+ = -0 . I. 2 7 4 + J ) 

CONSTANT C 105=0 . 6 9 7 ( C- C 2 , C i C6=- 0 . 7 2 C l 0 -C 3 ,C LC 5 = -0. ?4 3 3F- Cc 

G - OIN S T AN T C 1 0 9 = 0.2 l v E-- . Q I 

CONSTANT DU=- 3. SAO^t-C L ,!) 5 ? = C. ; A j'F+T 1 , ^3 3 = 3 . t ,<:7 0 - 3 1 
CONST ANT D44= -C . 3 7 50 c- 0 3 , C "4=-0 . FO 0 l.? + C2,0(r. V = -0 . 77 1 

CUNoT ANT 0 72 - J . 7 , > 5 ^ "H J ,1. . C 7 > - . 1 7 i- 1 h J 1 T 4 "="T1 • . n ~ - ■+ ^ 

CON ST AN T U54=- 0 . 4 1. 31 6- 0 l , 13 1 0 3 = 3 . 7 2 5 I " - C ! 

CONSTANT -£1=3.0, C2=C .0 , ' C 3 =0 .0 , ~ 0 4=3 .C , ~C 5 = 3.0 

CuN S I AN I 3 C S = J . . t '• U 7=2 .0 > I L o = 2 . > ! . 5 = ' . C . 1 ( l v- ). i 

CONSTANT C X V = 0 .0 5, U>jA = 0 .0 ,07.0=0 .0 , V.: * =0.0 

OY NAMiC 

OR HO 5 C = C U * n 4 H U S + C 1 2 '•= 0 T S + C l 3 * ^ p R + O ). 1. -h r X v 
ORHOS=IMTGkL ( - C !. , I3R hCS C ) 

b T S i ) -c 2 1 ’ HU K NO 5 itrg-j r wr Hrt e±5 * E"T 5 N-Hrft 3 + r P * 

n T S = 3 NTGRK “ C 2 , C T s 0 ) 

OT S W 0 =C 3 I * D A HO S + C 3 ?. * ; )T S + C 3 3 « 1 T S W t-r - h i- ” o p +q 3 7 * 0 0 Hi A 

D 1 SU=IN I GKirn C 3 , n I S 1 -'.) ) 

UX0=C4iv0kH:)S +C42v-UT S +C 44 + 0 X +C 45 +')'•» +0 4 6+CU‘.?+C A7*0 TP'J . . . 

— I i-L 4.5 It DP j-i tC-4-a ;i- N T ift.t: 3 4 A 2- 10 S : 

nx = I NTGRL ( f C4 .0X0 

0 * 0 = C 5 ). * l) k h 1 ) S + C 5 2 1 J ' T S + C 5 4 * r ' X + 0 5 5 * ' 1 1 + C 5 £ <3 w *• C 5 7 • = n TO'*' ... 

1 1. ? r v tn * — 1 C8*7 ‘ . ‘ w T/ « 1 i 5 r r*i*{7~H — 

n,i = I NTGRL ( JC5. CM/ ) 

0'.* : /jo = C6l*DRHUS+C62vorS+C6 4*nx + C35- , ‘O/. , + 0 0( « hCc7^0To -I ... 

+(. !),v:, in t. , ‘ : 

l)NW= i TiT GKL ( i'C6 , CC VJ) 

0 t a k n =c 2 3 a- n 1 s ± r 7 7 f 1 a < , k + 0 2 a ± c r- * + 0 1 2 a r *, r + r 7 2 j - "na 

UTR w= TNTGk i. ( - c 7 , 0 r r, v.rj ) 

npriD=C3l=OPHOS + C8 2»*‘i TS +C6 4*C X+O.TT •' ' + 1 0 f N ) 1 - 1 » h f. 8 4 o J h C 3 F ; * ; T h ... 

t-Ht-d A^Ow 

00 8= 1 NT OP !. ( 0>; , i.p H' ) 

OT v/0 = C9 -*■ * DX + CO 5 * 0 -'i + C -6 h 0.5° *00'’ +0 94 >)r T Wf v-.' Ay-nc 

'll Ai= N TTF L ( : C C , T M* ) ~ 

UY 0= C 10 + «CX +c to 5 )•• HC lot A U40 1 no '-op A f r 1 r q-j T ,. + n 1 .14-AO . 

. 3 V = T N TO , . 3 I ( : C l a . r Y , 3 1 

OP S= A o 3 ->0 T S + A 1 j 1*0- or s 

DU 8 =n 73 ■* OP B -r> 7 <* <■ TT S- 87 5 * 0 K H S 

*■ ) ' . * 1 5 = 11 7 6»‘ CX Y * ' 3 7 7 ■' 1 . 1 >; i-’ -. - rd r > '-*• r- 6*8 



103 



PAGE HEj:GHT-3 , !» ’■ 071-1=5 
page xyplut 



-QO r PUT T l V h > S 

L A d fc i_ OP HO S V A T ' ' 

output T'lMr-.ms 



Abet. OTS VS T . V |- 
OUTPUT TJME.OTSW 



-u abfl ax-s 



OUTPUT 
LABEL OX 



— V > 

T T E E « 0 X 



VS T 



" OUTPUT r>'B, 
LABEL Ow VS 
OUTPUT TYM E. 



ttt 

T M c 



f)v 

T 



>{ 

r777- 



L A ti 1 . 1 . U W V V S 
OUTPUT T < ME *0 TB ■/> 
- 1 ABEL — C4-&1 — VS T 



OUTPUT T.ri'-E t OPB 
LABEL UP B VS T m f 

-O utput T' 



OT7 OTu 
L A B E L 0 T L VS T " P r 



OUTPUT T 



UY 



LABEL OY VS T-WE 
OUTPUT TIME, OPS 
-L A B r l, — y-P-S — V-S 7 i 'j r 



OUTPUT T ■ y F , 0 i w 
LABEL Ow S VS T * (•-: 
- output 



LABEL OWB VST */ p 

T^MEP E ■ UT M = 300 .0 , OUT DEL =3 . 0 tPTOE' -T.C 

CONST ANT 0 W c = } . 0 3 _ « n X V = U . C » C k A - 7 . ),(-)'<•" - 3 . 
■ 



CONSTANT CW= 0 . 
EN D 

-7 ~ ; TUP 



iD B\=u . G » C ri r =0 .0 , OX V- 7 



104 



BIBLIOGRAPHY 



1. Chien, KL, Ergin, El, Ling C and Lee, A, " Dynamic 
Analysis of a Boiler", Trans A ma r So c Mech Enqr s, v. 
80, p. 1809- 1819, 1958. 

2. Whalley, B, The Control of M arin e Prop ul sion Plant, p. 

1.1-4.28, PhD Thesis, University of Manchester, 
1976. 



105 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Library, Code 0142 2 

Naval Postgraduate School 

Uonterey, California 93940 

2. Department Chairman, Code 69 1 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

3. Assoc. Prof. Thomas M. Houlihan, Code 69Hm 5 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

4. Assoc. Prof. Robert H. Nunn, Code 69Nn 1 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

5. LT Chusakdi Senanikrom, RTN 3 

Education Section Personal Department 

Royal Thai Navy, Bangkok Thailand 

6. Defense Documentation Center 2 

Cameron Station 

Alexandria, Virginia 22314 



106 






















T hes ! s 

S **1735 

c.l 



Senanikrom 



76931 



si!" s 






Thesis '» 7 1- q Q 1 

S41735 Senanikrom 

State variable anal- 
ys i s of a boiler sys- 
tem. 



\ 



lhesS41735 

State variable analysis of a boiler syst 




3 2768 001 94490 3 

DUDLEY KNOX LIBRARY 



